Local elasticity map and plasticity in a model Lennard-Jones glass. 
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In this work we calculate the local elastic moduli in a weakly polydisperse 2D Lennard-Jones 
glass undergoing a quasistatic shear deformation at zero temperature. The numerical method uses 
coarse grained microscopic expressions for the strain, displacement and stress fields. This method 
allows us to calculate the local elasticity tensor and to quantify the deviation from linear elasticity 
(local Hooke's law) at different coarse-graining scales. From the results a clear picture emerges 
of an amorphous material with strongly spatially heterogeneous elastic moduli that simultaneously 
satisfies Hooke's law at scales larger than a characteristic length scale of the order of five interatomic 
distances. At this scale the glass appears as a composite material composed of a rigid scaffolding and 
of soft zones. Only recently calculated in non homogeneous materials, the local elastic structure plays 
a crucial role in the elasto-plastic response of the amorphous material. For a small macroscopic shear 
strain the structures associated with the non-afhne displacement field appear directly related to the 
spatial structure of the elastic moduli. Moreover for a larger macroscopic shear strain we show that 
zones of low shear modulus concentrate most of the strain in form of plastic rearrangements. The 
spatio-temporal evolution of this local elasticity map and its connection with long term dynamical 
heterogeneity as well as with the plasticity in the material is quantified. The possibility to use this 
local parameter as a predictor of subsequent local plastic activity is also discussed. 

PACS numbers: 



I. INTRODUCTION 

It is commonly acknowledged that the mechanical 
properties and the rheology of a wide class of amor- 
phous glassy materials involve localized structural re- 
arrangements of the order of five interatomic distances 
involving of the order of 100 particles in three dimen- 
sions or 20 in two dimensions [l|, These events 
are frequently compared with the rearrangements that 
take place in aging glassy materials where local events 
are thermally activated [1, 0]- It has been proposed 
that these rearrangements can organise during a me- 
chanical deformation, through a cascade mechanism to 
form shear bands [H, d, H B Q- This phenomenon of 
strain localization has been observed experimentally in 
alloys, metall ic g lasses, polymers, granular media, foams 
and colloids [13 [O, [H as well as m numerous simu- 
lations, both of model systems such as Lennard-Jones 
glasses [E I3I [T^ . [isj . as in more realistic simulations 
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On a theoretical level various current 
models of the rheology of glasses predict reasonably well 
the macroscopic mechanical properties of these materi- 
als i, S [10, [in, [13, [H, [H, [li. Hi, [13, [H, [H, [M, [Ml, [13, 

[Hi [li I . These models generally involve the consideration 
of zones with prescribed elasto-plastic properties, whose 
microscopic identification remains however elusive. So 
despite a recent and rich literature concerning the con- 
nection between the structure of the glass, the intrinsic 
dynamics of these irreversible events and the dynamical 
heterogeneity in amorphous systems d, [13, [H 113; HI; 
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M IM IM l39|, |4d, |4l|, |4a IM |4J; |45| questions such as 
what type of, where and how local plastic rearrangement 
occur in a deformed glass remain unanswered. 

A common general idea is that the deformation will 
take place in 'weak' zones, somehow characterized by 
abnormally low elastic constants and increased mobil- 
ity. This picture of the glass as composed of a patchwork 
of 'rigid' relatively strongly bonded (but amorphous) do- 
mains separated by 'soft' regions (walls) has long been 
postulated ^fi] ^^id is at the heart of many theoretical 
models of the glass transition [13]. Moreover within this 
framework it appears tempting to relate the dynamical 
heterogeneities observed in glassy materials near and be- 
low the glass transition temperature to the spatially inho- 
mogeneous elastic constant network. Experimental evi- 
dence of the heterogeneous glass structure is shown in [ill 
and models [1^, |50| including fluctuating elastic moduli 
have been introduced in the last two decades to describe 
some acoustical (boson peak [sij ) and thermal (specific 
heat anomalies [l^ilHillJ) properties of glasses. In order 
to check such assumptions, and attempt to relate struc- 
ture and dynamics in amorphous materials a new route 
was recently opened by calculating local elastic constants 
in such heterogeneous material [55| . Simulation methods 
for the calculation of the elastic constants of solids can be 
classified in two groups. One can either use equilibrium 
fluctuation formulas, as described in [s^ or perform ex- 
plicit deformations and derive these constants from the 
stress-strain response of the material. In this paper we 
implement the deformation approach in a framework of 
continuum mechanics developed in , that extends the 
applicability of classical continuum mechanics to smaller 
scales. While this method requires several computations 
to calculate the different elastic constants, it remains sim- 
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pie and accurate, and has the advantage to give a sys- 
tematic estimate of the deviation from hnear elasticity as 
a function of the coarse graining scale. 

Section[TT]presents the system under study, a well char- 
acterized two dimensional Lennard-Jones system, and 
the method used to obtain the local elastic moduli. The 
results for the moduli and the resulting picture of the 
glass structure are described in section IIIII Finally, in 
section HVl we show how the local elastic inhomogeneities 
are connected to the plastic activity as our system un- 
dergoes quasistatic shear deformation, and the relation 
between local elastic moduli and local dynamics is also 
explored. 



II. METHODS 

A. Sample preparation 

The systems studied throughout this paper are slightly 
polydisperse two-dimensional Lennard-Jones glasses, de- 
scribed in detail in [57| . The interaction energy between 
two particles is 4e((2^)i2 - (^^)^) with cr, the di- 
ameter of particle i. Our typical sample contains 10 000 
particles in a square simulation box, with a polydisper- 
sity of 5% on the size of the particles. All numerical val- 
ues in the following are expressed in Lennard-Jones units 
(LJU), where the average particle diameter and the in- 
teraction energy are equal to unity. The corresponding 
density of the glass p — 0.925LJU was chosen to mini- 
mize the initial pressure of the quenched amorphous solid 
to P = 0.25. The quenching procedure, from the liquid 
state to well below the glass transition temperature, that 
produces the initial configurations consists in a sequence 
of thermal steps - at temperatures T = 1.0, 0.5, 0.1, 0.05, 
0.01, 0.005 and 0.001 - each lasting 100 LJU and ther- 
mostated through simple velocity rescaling. This quench 
is followed by an energy minimization scheme to bring 
the system at a local energy minimum and at zero tem- 
perature, this protocol is also described in detail in [57| . 
Two protocols were used to shear the material. In pro- 
tocol one, two layers of particles with thickness 2 LJU 
are singled out and assumed to constitute parallel solid 
walls that will impose the deformation to the system. 
The resulting cell has a thickness Ly = 100 and a width 
Lx — 104. The configuration is submitted to shear by 
applying constant displacement steps Sux to the parti- 
cles of the upper wall, corresponding to an elementary 
strain of about e = 5.10~^. Between two displacement 
steps the entire system is relaxed with fixed walls into 
its nearest equilibrium position. The total strain of the 
sample under these rigid walls boundary conditions is 
itotai = 1-65 (i.e. 165% deformation). In protocol two 
the system of thickness Ly = 104 and width Lx =104 
is sheared with Lees-Edwards periodic boundary condi- 
tions. The elementary strain step during the shear is here 
e = 2.510~^. The total strain under these Lees-Edwards 
boundary conditions is etotai = 0.5 (i.e. 50% deforma- 



tion). To check possible size effects a series of glasses of 
different sizes (containing up to 250 000 particles) were 
produced under the same quenching procedure and ana- 
lyzed in the very early linear domain {ttotai < 10~^). 

In the following section we derive the coarse- grained 
(CG) method used to calculate the local elastic tensor C 
as well as other local CG fields such as the strain tensor e, 
the stress tensor cr, the density p. This method is applied 
to configurations in the transient regime e < 2.510^^ as 
well as in the fully developed plastic regime e > 2.510^^. 
The domains analyzed are shown in figure [TJ 
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FIG. 1: Stress-strain mechanical response of the 2D polydis- 
perse Lennard-Jones model glass under shear, indicating the 
separation between an elastic loading phase and a plastic flow 
behavior. Rigid (protocol one) and Lees Edwards (protocol 
two) boundary conditions are respectively represented in thick 
and thin black lines. 



B. Measuring local elastic constants 

As stated in the introduction, various methods have 
been proposed in the literature for the numerical cal- 
culation of elastic constants. The elastic constants are 
defined [111 as the coefficients Cap-yS of the second order 
expansion of the strain-energy density as a function of 
the local Green-Lagrange strain components 

W-Wo_^,o 1^ 

y — L^Q/jEa/S + 2^a/375^a/3e75 + (1) 

or as the first order expansion of the local stress compo- 
nents (Tap as a function of the local linear strain [soj 

o'a/3 = C^p + Cafjjse''"s -f ... (2) 

with 

1 dW 



This second definition can be seen as a linear local fit of 
the energy and is known as Hooke's law. In case of an 
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unstressed solid, the two definitions are strictly equiva- 
lent, but in case of a solid with initial stresses {C° ^ 
as in our Lennard- Jones glasses), the difference between 
the two components Cap-yS and C'^p^g depends on the 
quenched stress components, due to the non-linear de- 
pendence of e as a function of e''" [s^. In our case, we 
have checked that the contribution of quenched stresses 
can indeed be neglected [B^I, but this formal difference 
must be mentioned. 

At finite temperature one can use a fluctuation type 
method to determine ([1]), and for a system at equilibrium 
the elastic modulus tensor and the stress fluctuations are 
related through the following fluctuation formula [6l|, [g^. 



a 



f^Born 
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[{iapt^s) — {tal3){t-ts)\ ■ 



(3) 



The term in square brackets (here the brackets denote 
a thermal average) is called the fluctuation term of the 
microscopic stress tensor t^g- The Born term Cj^^^g cor- 
responds to the instantaneous elastic modulus for a sam- 
ple under an uniform and infinitesimal strain. Equation^ 
has been used for a long time to compute elastic constants 
of amorphous materials [64] but only recently applied to 
the calculation of local elastic constants in amorphous 
polymers ^] , in metallic glasses [g^ , in composites 
as well as in small samples of Lennard- Jones glasses 
above the glass transition temperature. If one is inter- 
ested, on the other hand, in the mechanical properties 
of a material at zero temperature one must use the zero 
temperature limit of equation [3l that was shown in [6^ 
to be 
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can be understood as the forces which would result from 
an elementary homogeneous deformation of all the par- 
ticles in the strain direction rj^s and where H is the dy- 
namical matrix of second derivatives of the potential en- 
ergy with respect to atomic positions. At zero temper- 
ature, the 'relaxation-fluctuation' term (second term on 
the right hand side of equation [4]) was shown [s^] to ac- 
count for an important fraction of the absolute value of 
elastic constants in amorphous systems (this is also the 
case in crystals with a complex unit cell) . This failure of 
the Born term alone (first term on the right hand side of 
equation [3]) to accurately describe the mechanical prop- 
erties of the material can be traced back to the existence 
of a non-affine deformation field, which stores part of the 
elastic energy. Unfortunately, the direct evaluation of the 
'relaxation-fluctuation' term is not straightforward as it 
necessitates the inversion of the Hessian matrix. Hence 
the actual procedure to accurately determine the elastic 



constants of athermal systems generally consists in car- 
rying out explicitly an affine deformation of all coordi- 
nates and a corresponding deformation of the simulation 
cell, then letting the atomic positions relax towards the 
nearest energy minimum within the deformed cell. One 
calculates the corresponding local stress increments and 
linear strain components, and can subsequently recover 
the elastic tensor coefficients by solving Hooke's law [51 

In this paper we propose to extend this method 
to the calculation of the local elastic constants at 
different scales of coarse graining. Following [5^, we 
present first the coarse-grained expressions that we 
use to measure locally the stress and strain fields, we 
then pursue with the derivation of the local elastic tensor. 

Local strain tensor 

The simplest and crudest approach is the so called 
Voigt assumption where one assumes that the strain is 
homogeneous in the sample. While this approach is ex- 
act for simple unit cell crystals it cannot be valid for 
disordered materials and leads to wrong estimates of the 
elastic constants as discussed above. It was proposed in 
refs [20| as an improvement on the mean field assumption, 
that the local strain should be determined from a best fit 
of the actual particle displacements within a small sub- 
system to those generated by an adjustable strain tensor. 
Other possibilities for defining a local strain include the 
consideration of a dyadic tensor built on the links be- 
tween neighboring particles ,£8,]. In this study we will 
use an alternative definition proposed in according 
to which the (linear) strain tensor is written in terms of 
a coarse grained displacement field: 



lin 



9<"(r,t) 



(6) 



where the superscript 'lin' denotes the linear order in 
the displacements. The linear displacement u''"(r, i) is 
defined as the linear dependence on the displacement of 
the individual particles, of the temporal integration of 
the coarse-grained velocity. The latest is computed by 
the way of the momentum density whose coarse-grained 
expression satisfies the mass conservation equation. As 
shown in ref [56j, by integration of the coarse-grained 
velocity field u(r,t), one obtains the following expression 
for the linear order in displacement: 



lin 



(r,t) 



X;,-mj(/)[r-rj(t)] 



(7) 



where 0(r) is the coarse-grained function, for example 
as in all calculations done in this paper, a gaussian of 
width W. u''"(r, i) is a continuous and differentiable 
function of space that allows to compute e''" by spatial 
derivation. The difference u''"(ri,i) — Ui{t), where 
u,; is the actual displacement of the particle i, is the 
fiuctuating part of the displacement [69|. It can not 
be negligible when the displacement of the particles is 
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strongly varying with space, as it is shown at very small 
scales in our Lennard- Jones glasses [69j . 

Local stress tensor 

The most commonly used definition is the atomic stress 
introduced by Irving Kirkwood 70] as 



and one has in 2D the relation T = CE that means: 



(8) 



Vi being the volume of the Voronoi cell associated with 
atom i, and fi/ja the force exerted by atom j on atom i 
in the direction a. If the velocity term is neglected this 
formula when averaged over a sufficiently large sample 
corresponds to the the Cauchy stress (which measures 
the actual mechanical force per unit area), at smaller 
length scales though it does not, strictly speaking, verify 
the momentum conservation equation. 

An expression for the stress components that satisfies 
the momentum conservation in the framework of the pre- 
vious coarse-graining has been obtained in [56| 



a„^(r,i)= (9) 
~ ^ 51 XI -^Vj" i*y^]P (*) /o ^^'t> ~ (^) + (0] • 

with (j) the coarse-graining function. This expression has 
also been used recently in [7l| . 

Derivation of linear elasticity - local elastic 
tensor 

Now having calculated the stress and strain tensors 
locally we are in a position to derive the corresponding 
components of the local elastic tensor as in the case of a 
macroscopic deformation. This presuposes that equation 
[2] remains valid at a local level and one can therefore ap- 
ply the same symetry arguments as in the macroscopic 
case to reduce the number of independent coefficients of 
the elasticity tensor. Since the linear strain is symetric by 
definition, and the stress is symmetric in the absence of 
torques there are at most 9 constants in 2D. As discussed 
in [5^, [5^ the existence of a strain energy function from 
which the equations of elasticity can be derived by vari- 
ational methods implies a further symetry of the elastic 
tensor : Cap-yS = C^sap , reducing the number of indepen- 
dent thermodynamic stiffness to 6 in 2D. To extract the 
6 independent elastic coefficients necessitates at least two 
independent deformation modes on our sample. Each de- 
formation provides three linear equations for the moduli. 
The general stress-strain relation in terms of matrices is 
written as follow, using a Voigt type notation. For each 
coarse- graining scales W , the coarse-grained stress and 
strain components are measured on a grid, and are ex- 
pressed respectively as the 3*1 column vectors T and E 
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This expression can be compared with the expression ob- 
tained in the framework of homogeneous and isotropic 
linear elasticity: 



^xx 

^yy 
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■ (11) 

where fi is the shear modulus and A the Lame coefficient 
((A + /i) is the inverse of the compressibility modulus in 
2D). 

The use of 2 such deformations therefore closes the 
system of unknowns, giving 6 equations for 6 unknowns. 
Nevertheless, in order to estimate the deviation from lin- 
ear elasticity the numerical procedure used here consists 
in applying three different uniform deformation modes, 
two uniaxial stretching parallel to the x and y axis and a 
simple shear parallel to the x axis. This procedure pro- 
vides 9 linear equations for the moduli. The stress com- 
ponents which are not used in this procedure are then 
calculated using these elastic moduli, and their values 
compared to those computed directly using equation [9l 
As a measure of the extent to which the system is de- 
scribed by linear elasticity at a given position and for a 
given value of the coarse-graining scale W we use the root 
mean square. A, of the relative differences between the 
stress components calculated by employing the measured 
moduli and the directly measured exact values (normal- 
ized by the norm of the exact values). 

For each configuration we calculate also the 3 eigen- 
values Ci and eigenvectors Ei for i — 1,2,3 of the local 
tensor C. The comparison between eg nations 1101 and 1111 
would provide ci — 2^, C2 = 2/j, and C3 = 2(A + //) in an 
homogeneous and isotropic system. We will now discuss 
the results obtained in our model Lennard- Jones glass, 
as a function of the coarse-graining scale W . 



III. ANALYSIS OF THE LOCAL MODULI 

The domain of validity of Hookes's law can now be 
measured by the H^-dependence of the error function 
A. It is shown in figure [21 This figure shows that 
the error A goes progressively to zero, thus validating 
Hooke's law at large coarse-graining scales W . However, 
this convergence obeys a power-law A c>c W'"^'^^ and 
therefore doesn't exhibit a characteristic scale above 
which Hooke's law would be true. Nevertheless, we can 
see in figure [2] that the error A is already less than 1% 
for W >'b. It means that above = 5, the error made 
in computing the stress components using the elastic 
modulus is less that a hundred of the actual value of 
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FIG. 2: Deviation A from linear elasticity as a function of 
the coarse-graining parameter W. For W > 5 Hooke's law is 
satisfied locally with more than 1% accuracy. 
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FIG. 3: The three eigenvalues of the elasticity tensor averaged 
spatially as a function of the coarse-grain length W. Inset: 
log-log plot of the convergence to the limit value obtained by 
a coarse graining on the whole system size. 
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FIG. 4: Anisotropy parameter na as a function of the coarse- 
graining scale W . 



the stress, that is far smaller than the contribution of 
quenched stresses for example. We can thus consider 
that the system obeys Hooke's law reasonably well on 
scales larger than W = 5. Below that scale, different 
factors could explain the departure from Hooke's law: 
the first one is the coarse-graining. Below W — 5 
it has been shown in a previous paper [6^ that the 
contribution of the coarse-graining deformation to the 
actual one is small. The contribution of the fluctuating 
field can not be neglected, giving rise to high values of 
the real strain and therefore a significant decrease of 
the elastic moduli. In fact, the fluctuating field is not 
differentiable, so that it is not possible in this case to 
compute quantitatively the linear strain components. 
It is one of the interests of the use of a coarse-graining 
field to deal with differentiable displacements fields. 
The fluctuating field appears thus as a 'noise field'. 
An additional strongly fluctuating term should then be 
taken into account for < 5 in an attempt to describe 
accurately the mechanical behavior of the material. An- 
other contribution to the departure from Hookes's law 
at small scale is due to the coupling to second and third 
neighbors outside the volume element. This contribution 
(not taken into account here since we are restricted to 
first order derivatives in the displacement field) could 
be introduced by considering the contribution of higher 
order derivatives of the deformation in the framework 
of linear but long-range elasticity [T^I- This paper is 
devoted to the measurement of elastic moduli, thus we 
will leave these considerations for further studies, and 
focus now on the computation of elastic moduli in the 
domain of validity, that is for > 5. 

Figure [3] shows the average value of each of the eigen- 
values < Ci > of C* as a function of the coarse-graining 
scale W. The notation c stands for a spatial averaging 
over the sample and angular brackets stand for a statis- 
tical averaging over different configurations. The spatial 
averaging is obtained by computing the elastic moduli on 
a grid with elements of width W/2. The number of the 
different samples used in the statistical averaging - for 
each given size L - is the same as indicated in the figure 
[21 We get flrst the average value of the (L/(VF/2))^ values 
obtained on the grid, for each sample. Then we average 
the values obtained on the different samples. The flgure[3] 
shows a progressive convergence for large W to the values 
obtained in the framework of homogeneous and isotropic 
linear elasticity. Indeed, for very large W, < cT > and 
< C2 > go to the same value 2/i « 22 obtained also by 
looking at the macroscopic response of the system to var- 
ious mechanical soUicitations, and < 7^ > converges to 
twice the inverse compressibility 2(A-|-/z) « 102 measured 
as well by the global response of the sample [13, [ill . The 
method used here to compute the elastic moduli of the 
system for large W appears thus to be consistent with 
measurements of the global response of the system. 

The convergence of the average local elastic moduli 
to their macroscopic value is independent of the sys- 
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tern lateral size L as long as, W < 0.5L. For larger 
W ^ the boundary conditions (Lees Edwards or fixed 
walls) may affect the convergence and cause finite 
size effects. As shown in the inset of figure [31 the 
moduli decay approximately as to their limit 

value, with a ~ 0.87 for ci and C2, and a ~ 2 for C3. 
The inverse compressibility converges thus far more 
quickly to its macroscopic value, than the shear moduli. 
The difference between the macroscopic value and the 
spatial average of the small scale measurements of 

« 22 while < cr(W = 5) >w 18 for example) is 
due to the inhomogeneous strain field. By looking at 
(< C > — Coo)/Coo (inset of figure [3]), it appears that 
for > 5, the discrepancy to the macroscopic value is 
already less than 1% for < C3 >, while it becomes less 
than 10% only for W > 10, for < ci > and < C2 >. 
We conclude that for 5 < < 10 the moduli are well 
defined, but the measured values are not compatible 
with homogeneous elasticity since the different moduli 
involved at different coarse-graining scales have different 
scale dependence. We did not find any solid explanation 
for the non-trivial power-law appearing in this conver- 
gence. It appears that the convergence of the inverse 
compressibility is inversely proportional to the volume 
W^, and the convergence of the shear moduli closer to 
a surface effect cx W. Note also that one of the shear 
moduli is smaller than the limiting value, while the 
other is larger. This difference between smooth and hard 
directions will now allow us to quantify the anisotropy 
of the local mechanical response. 

The anisotropy measured at small scale can be quan- 
tified by the ratio (< ci > — < cT >)/(2^) that goes to 
zero for large W. We call it the anisotropy parameter iTa- 
It is shown in figurelH It can be noticed that ci and C2 al- 
ways obey < ci >< 2^ and < C2 >> 2/j,, so Wa > 0. The 
decay of the anisotropy parameter iTa (figured]) obeys also 
a power-law oc l/W^-^-^, close to 1/W. These power-law 
decays prevent us to define properly a characteristic scale 
above which the homogeneous and isotropic behaviour is 
recovered. In figure UJ we see that the anisotropy param- 
eter becomes less than 10% for W > 20 only. It means 
that, below W = 20, it is possible to find locally a well 
defined direction associated with a very low local shear 
modulus. At this scale, the anisotropy in the mechanical 
response cannot be neglected. 

The preferred direction of strain is given locally by the 
analysis of the eigenvectors Ei , E2 and . Each eigen- 
vector contains the 3 distinct elements of a 2D strain 
tensor whose eigendirections ei and 62 (ei and 62 arc 
orthogonal) are computed. We plot on figure [D the dis- 
tribution of the local quantity Si = {tr{Ei))^ /2tr{Ef), 
that takes the value if the deformation is pure shear and 
1 for pure dilatation. One observes as expected that the 
two deformations associated with the two lowest eigen- 
values are of pure shear type while the third deformation 
is a pure compression. 

In order to explore more deeply the inhomogeneities of 
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FIG. 5: Distribution for each eigenvalue Ei of the local elas- 
ticity tensor of the quantity Si = {tr{Ei)f /2tr{Ef). The 
distribution is obtained over the entire plastic flow regime. 



the elastic moduli inside the system, we will now study 
their distribution as a function of W, and then their spa- 
tial correlations. The distributions of ci, C2 and C3 are 
shown in figure [6] for various W and for N = 216225. 
First, we can see in these distributions, that zones with 
negative moduli can appear if the coarse-graining scale 
W is sufficiently small, as already observed in [5^. It 
is not in contradiction with the mechanical stability of 
our system, since the local elastic moduli computed here 
are only part of the second order derivative of the to- 
tal mechanical energy, due to the coarse-graining, as well 
as to the non trivial dependence of the non-affine local 
deformation as a function of the applied displacement. 
The rescaling of the distributions by W is also shown 
on the figures [HI It is very good for sufficiently large 
values of W (typically W > 10). The variance of the 
distribution as a function of W is shown in the inset of 
figure [BJi. It decays as 1/W and is always smaller than 
the corresponding average value. It is thus impossible to 
identify a characteristic lengthscale by the comparison of 
the variance and the average value of the moduli Ci . 

The decay of the relative fiuctuations Ac/ < c >oc 
1 /W for a given W can be interpreted in the framework of 
a sum of uncorrelated variables with finite variance, the 
distribution being nearly gaussian. The apparent rescal- 
ing of the distribution thus corresponds to a sum of spa- 
tially uncorrelated variables. Note that, while this ratio 
is very small for C3 for every value of W (AC3/C3 < 0.01 
for > 5), it is much larger for ci (Aci/ci > 10% while 
ly < 15 for iV = 216225). We can thus conclude that 
the inhomogeneity is far more pronounced in the shear 
modulus ci than in the compressibility C3. The inhomo- 
geneity of ci-modulus is even far from been negligible 
while W < 15. 

By comparing this result with the result obtained for 
the isotropy, we can conclude that at a scale W > 20 the 
system becomes reasonably isotropic and homogeneous. 
Below this scale, it is homogeneous but not yet isotropic 
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FIG. 6; The figures (a), (b) and (c) show, respectively, the 
distributions ol the three lowest eigenvalues ci , C2 and C3 for 
different coarse-graining lengths W . The inset of each figure 
show the distributions rescaled by W. Figure (d) shows that 
the distributions get narrower as W increases to converge to 
the macroscopic values of the elastic constants. The variances 
of the distributions is plotted as a function of W . The vari- 
ances normalized by the average values of the distributions 
are also plotted. Both curves show a power law evolution 
with W. 



for 15 < W < 20. All these results are summarized in 
Table m Of course, these values must be nuanced by the 
fact that the criteria used are a little arbitrary and sim- 
ply related with a comparison of orders of magnitude. No 
characteristic lengthscale can be clearly identified since 
all the quantities checked here have a power-law depen- 
dence with the scaling length W. 



w 





5 


10 


15 


20 


Hooke's law 


NO 


YES 


YES 


YES 


YES 


Homogeneity 

2^ <io% 


NO 


NO 


YES 


YES 


YES 


Ac 

7=r< 10% 

(c) 


NO 


NO 


NO 


YES 


YES 


Isotropy 

2^ < 10% 


NO 


NO 


NO 


NO 


YES 



TABLE I: Analysis at different coarse-graining length scales 
W. 



The distribution of the elastic moduli has been checked 
along the full deformation process. We can thus compare 
the distribution of moduli during elastic deformation and 
during plastic flow. We can also compare these distribu- 
tions before and after a plastic rearrangement occurred 
in the system. The plastic rearrangements are identified 
here as in [T, ^T^l by the decrease of the total stress as a 
function of the applied strain. We see in figure[7]that the 
distribution is progressively displaced to smaller values of 
the shear modulus, before it reaches a plastic plateau, but 
remains globally unchanged during all the plastic flow. 
During the plastic flow, the difference between the dis- 
tributions appears on extremal values. Before a plastic 
rearrangement occurs, the smallest value of ci (open cir- 
cle in figure [8|) is generally smaller than after the event 
occurred (full circle). The distribution of the smallest 
value of this elastic modulus is shown in the inset of fig- 
ure [5] for the full deformation process. It confirms that 
the smallest values are smaller before a plastic rearrange- 
ment than after. We will come back to this observation 
in section HVl 

Finally, it can be interesting to analyse in details the 
spatial correlations of the moduli in our systems. The 
spatial correlations of the lowest modulus ci are shown 
in figure [5] for various W and N — 216225. The spa- 
tial correlations go to zero at large distances. It shows 
spatial oscillations with very small amplitude, that are 
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visible while W < 10, but that disappear for W > 
10. The distance between successive maxima is about 
a few tens interatomic distances, but seems to be size- 
dependent. Unfortunately, our data are not sufficiently 
precise to allow us to characterize this size-dependence. 
For r < 3W only, the spatial correlations arc con- 
trolled by the W^-dependence of the coarse-graining func- 
tion (see figure [IJb)). It can be fitted by a gaussian 
cx exp{—{r/W)'^/1.7). We can thus conclude that the 
spatial correlations are dominated by the W^-dependence 
of the coarse-graining function at small distances, but 
displays oscillations at larger scales. The domain for 
which these oscillations are visible {W < 10) corresponds 
to the domain in which the heterogeneity in the distribu- 
tion of the moduli is noticeable. All these results show 
that a coarse-graining at scales W >> 10 will loose in- 
formations (on the heterogeneities, on the spatial corre- 
lations, even on the local anisotropy.). In the future we 
will thus use only the value W = 5 for the coarse-graining 
scale. 
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FIG. 7: Distributions of the local shear modulus ci averaged 
over the elastic regime and for two configurations in the plas- 
tic regime, just before and just after a rearrangement. In- 
set: Evolution with the macroscopic strain of the instanta- 
neous macroscopic shear modulus defined as Aaxy.macro = 

1-imacro A^xy , macro • 



IV. STRUCTURE AND DYNAMICS 

A. Structural relaxation 

In the previous section we have shown how the sys- 
tem could be decomposed into regions of different elastic 
stiffness. We now discuss how this elastic heterogeneity 
is related with the 'dynamics' of the system undergoing 
quasistatic, plastic shear deformation. To this end, we 
obtained the local elastic parameter ci for configurations 
of the sheared system separated by an incremental strain 
of Ae ~ 5.10"^, during a set of intervals each within a 




Q— Q Min(Cj) before plast. event 
Q— □ Max(C|) before plast. event 
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■ — ■ Max(Cl) after plast. event 
« — « Min(Cj) all steps 



FIG. 8: Evolution with the strain of the maximal and minimal 
values of ci. Inset; Distribution of the maximal and minimal 
values of the local elastic modulus ci over the entire plastic 
regime for all configurations, and for configurations just before 
the occurence of a plastic event. 




FIG. 9: (a) Spatial correlation function of the local shear 
modulus ci: ((ci(0) — cT){ci{r) — cT)). (b) Scaling of the 
function with the coarse-graining length W for small values 
of r. 



total strain of ~ 10% in protocol one, under rigid bound- 
ary conditions and ~ 6% in protocol two, under Lees 
Edwards boundary conditions. The coarse-graining scale 
W = 5 was chosen in all the subsequent analysis as the 
limit of applicability of linear elasticity. We recall that 
at this scale, < cT > (M^ = 5) w 18. 

One can first quantify the global relaxation time 
(strain) associated with the field ci by calculating the 
spatially averaged two point correlation function C(Ae) — 
(ci(r, e -I- Ae)ci(r, e)) where the notation A stands for 
a spatial average over the sample and angular brackets 
stand for a statistical average over the strain origins e. In 
order to relate the relaxation strain to the local rigidity 
of the material we also calculate the two-point autocor- 
relation function of the shear modulus conditioned by its 
initial value. For each sampled rigidity ci we plot the 
rescaled autocorrelation function C(Ae)/C(0)^/^. In fig- 
ure [TO] we see that the functions C(Ae)/C(0)^/^ tends in 
the limit of large strains to the limiting value cT, inde- 
pendently of the initial value of the shear modulus ci. 
Therefore in the stationary plastic flow regime the lo- 
cal shear modulus ci(r, £ + Ae) becomes uncorrelated for 
sufficiently large strains Ae from its value at the origin 
ci (r, e) , showing that in this model glass the local elastic- 
ity map does not phase separate into permanently rigid 
and soft regions but rather evolves dynamically under 
shear. In the inset of figure [TO] we associate with each 
sampled rigidity ci a relaxation strain treiax{ci) defined 
as the strain for which the rescaled autocorrelation func- 
tion has decayed by half. We see that ereiax(ci) is a 
monotonic increasing function of the local rigidity pa- 
rameter, that saturates for C(O)"'^/^ > cT- We see also 
that softer regions 

(C (0)1/2 < ^) relaxes more quickly 

than rigid ones. 

To our knowledge this result represents one of the 
first numerical evidence of a clear relation between a 
structural order parameter and the dynamics in a glassy 
system. The measurement of the local elasticity map 
presents the advantage to be independent of the speci- 
ficity of the glass under study, requiring only the mea- 
sure of a local stress and strain. This result confirms the 
description introduced in section IIIII of the material as 
a composite material made of 'soft' fast relaxing zones 
(for ci < ci) and of 'rigid' stable zones (for ci > ci). 
As seen in figure [TO] the strain associated with the rigid 
'scaffolding' of the material is found to be of the order of 
treiax — 1-5%. This valuc is similar to the strain neces- 
sary to enter the fully plastic regime tpiasUc — 2% (see 
figure [T]). As suggested in figure [TT] where the relative 
number of soft zones evolves in parallel with the total 
shear stress and reaches a maximum percentage (~ 60%) 
before a large plastic event, one can also see this typical 
strain as the necessary strain required to achieve perco- 
lation through the material of soft zones, i.e. when the 
material's rigid scaffolding is no longer connected p^ . 
Figure [TO] shows the evolution of the rigid scaffolding over 
a total strain of ~ 0.9%. 

This relaxation strain of the order of 1.5% can be com- 
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Autocorrelation function C(Ae) 



(ci (r, e -I- Ae)ci (r, e)) of the order parameter conditionned by 
the value of Ci(r, e) at the origin. The curves correspond from 
bottom to top to increasing values of ci (r, e) at the origin 
taken in the ranges [5-7,5], [7.5, 10], [10, 12. 5], 12. 5, 15], [20- 
22.5],[22.5,25],[25,27.5],and [27.5,30] as marked by the circles. 
The dashed line corresponds to cT. Inset: the typical relax- 
ation strain ereiax associated with each rigidity ci ~ C(0)^^^ 
defined as the strain for which the rescaled autocorrelation 
function has decayed by half. 




FIG. 11: Evolution of the total shear stress (thick line) and 
the relative number of soft zones (relative number of points r 
where ci (r) < cT) as a function of the applied shear strain, in 
th plastic flow regime. 



pared with the typical strain separating two irreversible 
rearrangements [l| in the sample Ace^ent ~ 0.1%. An es- 
timate of the number of plastic rearrangements required 
to renew the rigid 'scaffolding' of the material can there- 
fore be given as ereiax / ^^event — 15, hence typically 15 
events for 10000 particles. 



B. Relation between the local elasticity map, the 
local mobility and the long time dynamical 
heterogeneity 

We have characterized the dynamics of the underlying 
structure in terms of the shear rigidity order parameter. 
We want now to see how this structure is coupled to the 
displacement field in the sheared material. 

In order to describe the relation between the local 
elasticity of the material and its dynamics, one can try 
to quantify the connection between the local domains 
presenting small local modulus and an increased mobil- 
ity of the particles in these domains. A vast literature 
has grown in the last five years on the connection be- 
tween static structural properties and dynamical hetero- 
geneities in glass formers. While different approaches 
have partially failed to support such a link between lo- 
cal structure (such as local free volume, local inherent 
state potential energy, defects, Voronoi tessellation, local 
stress or strain...) and dynamics, Harrowell et al |42l | have 
recently shown that the spatially heterogeneous 'local 
Debye- Waller (DW) factor' (defined as the mean-squared 
vibration amplitude of a molecule over a time of approx- 
imately 10 periods of oscillation of this molecule) in a 
two-dimensional glass forming mixture could be mapped 
perfectly on the locally measured dynamical propensity 
that relates to the long term d yna mical heterogeneities 
in the material. Berthier et al (75| pursued this discus- 
sion showing how the influence of structure on dynamics 
is much stronger on large length scales than on shorter 
ones, and that the choice of the coarse-graining scale in 
the structure-dynamics problem is crucial. Here we make 
connection with this literature and claim that the local 
order parameter ci is a good candidate to establish a 
relation between structure and dynamics. One can un- 
derstand this assertion by the fact that ci and the DW 
factor contain a similar physical information in probing 
the local stiffness of the material. Of course in a qua- 
sistatic deformation one can not simply measure a local 
DW factor on short time scales and the order parameter 
ci is a good measure of the stiffness of a region. Follow- 
ing Cooper and Harrowell in [i^ we define a quasistatic 
analogue of the dynamic propensity as ([ri(e) — ri(0)] ) 
where ri{e) — ri(0) is the displacement without the afHne 
contribution due the macroscopic strain e. Unlike the 
original definition of the dynamical propensity in (43 | 
the average is taken here over all the particles in a given 
range of the order parameter ci and not over an iso- 
configurational ensemble of N-particle trajectories. Fig- 
ure [13 shows that the order parameter ci is indeed re- 
lated to the long term propensity and that soft regions 
(ci < ci) present an increased non-affine mobility in com- 
parison with more rigid regions of the material (ci > ci) . 

The relation between mobility (plasticity) and low 
shear modulus is illustrated in figure [14] by looking at 
how the spatial distribution of these two quantities are 
mapping onto each other. In this figure, the cumula- 
tive non-affine displacement field (that is essentially ir- 
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FIG. 12: Mean-square displacement (Ay^) for different 
values of the original local shear modulus ci. From 
top to bottom the curves correspond to ci taken in 
the ranges [0-2.5], [2.5,5], [5-7.5], [7.5,10], [10,12.5], 12.5,15], [20- 
22.5], [22.5,25], [25,27.5],and [27.5,30]. Inset: {Ay^} is shown 
as a function of ci for different relaxation strains. 



reversible) appears to nucleate from the initial reference 
state near the soft zones of the material and to grow in 
a cooperative manner up to the point where the mate- 
rial fails macroscopically forming a vertical shear band 
across the sample. In figure [T5| we distinguish between 
mobile and immobile particles (Fig. [T5| left), and be- 
tween soft and rigid zones (Fig. [T5l right). The mobile 
and frozen particles are identified somehow artificially by 
the amplitude of the transverse nonaffine displacement: 
Ay > 0.02 for a total strain of 1% and Ay < 0.02 re- 
spectively. The soft zones are identified by ci < cT and 
the rigid zones by ci > cf. In figure 1151 we plot the dis- 
tribution of shear modulus associated with each group 
of particles (mobile and frozen) and the distribution of 
nonaffine displacements for the rigid and soft portion of 
the sample. Figure [15] confirms the visual impression of 
figure [14] that most of the displacement is concentrated 
in the soft regions of the material and conversely that 
mobile particles are located in soft zones. 

We showed that most of the displacement occurs in soft 
zones. As seen in figure [HI the dynamics in the sheared 
glass is not trivial with regions that concentrate most of 
the non-affine displacement field and others that remain 
quiescent. The appearance of bursts of mobility seems 
therefore strongly dictated by the underlying heteroge- 
neous elastic structure of the material and one cannot 
understand cooperative dynamics in the glass without 
considering this underlying structure. We would like now 
to address the question of the degree of cooperativity of 
this mobility field in the material and its relation with 
the local elasticity map. In the literature the dynamical 
heterogeneity of ageing ([llj) or sheared ( [zll ) [z3| > Iz3 ) 
glassy systems is commonly quantified by a four-point 
correlation function defined as: 

X,{k,e)^^[{FAk,er)-{Fsik,e)f] (12) 
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FIG. 14: Wc represent the local map of the shear modulus at the same macroscopic strain values as in figure 13. These 
values correspond to the onset of plastic rearrangement of the material. These maps are superimposed with the non-afBne 
displacement accumulated from figure (a). Figures (a) to (f) correspond to a macroscopic strain of (a)2.5%, (b)2.55%, 
(c)2.65%,(d)2.7%,(e)3.25% and (f)3.4%. Note that the non-afBne field is multiplied by a factor 300 on fig (b) to illustrate 
the very strong correlation of the elastic non-affine field with the elasticity map for small incremental strain intervals (here 
0.05%). 
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FIG. 15; On the left the distribution of the local shear modu- 
lus for mobile and immobile particles. On the right distribu- 
tion of the nonaffine field during a total strain of 1% for rigid 
and soft regions. 



rigidity. One sees in the inset of figure [TBT a) that the re- 
laxation strain associated with each SISF grows hnearly 
with the local shear modulus below the average shear 
modulus (ci < cT « 18) and then reaches a plateau at 
a value of about Creiax ~ 0.85%. Again this provides a 
clear evidence of the connection between the structural 
order parameter ci and the dynamical response of the 
material. It cuts the sample into soft (ci < cf) and rigid 
(ci > ci) zones. 

In figure [TBFb) we see in X'liky = kp,t) that the 
number of particles evolving in a cooperative way grows 
linearly with strain from zero to ^ 500 particles for 
LE boundary conditions and to ~ 100 for rigid bound- 
ary conditions. It is indeed system size dependent and 
evolves as N.f{^/L), as shown in (goI. ItsI. [79l| . 

The maximum cooperativity is achieved at a strain of 
e ^ 2% in the LE case as well as in the rigid walls case. 
Theoretical predictions concerning the four-point corre- 
lation function are reported in ((41i]) where the authors 
focus their studies on static supercooled liquids near the 
glass transition temperature. Whether or not one can 
identify the dynamical correlation length scale to the 
typical spatial extent of the soft zones of the material 
remains unclear from this analysis and requires further 
studies. We note that the typical strain at which the co- 
operativity Xi{ky = kp,t) is maximal {emax ~ 2%) does 
not correspond exactly to the structural a-relaxation 
strain < 1% defined as Fs{ky — kp,ea) ~ 1/e cal- 
culated for the same wave vector, but the order of mag- 
nitude is actually the same for the systems studied here. 



C. Predicting plastic activity 



where (k, e) is the self intermediate scattering function 
defined by: 

F,(k, e)^J2 - "^^(0))) . (13) 

i 

It is important to note here that the symmetry of the me- 
chanical deformation introduces an anisotropy in k-space 
in the relaxation of the self intermediate scattering func- 
tion (SISF) Fs(k,e). Typically one has at the first peak 
(fcp) of the static structure factor a relaxation strain of 
Fs{kx = kp, e) of about 0.2% while for Fs{ky — kp, e) the 
relaxation strain is of about 1% (see figure frBFa)). This 
difference can be attributed to the formation of shear 
bands in our model glass, preferentially along the x axis, 
therefore increasing the mobility along this axis. These 
results seem in contradiction with recently reported sim- 
ilar studies in a model binary supercooled liquid, where 
an isotropic relaxation is reported for Fs (k, e) . It would 
be interesting to clarify, if these discrepancies could be 
attributed to the thermal agitation present in [t^I and 
absent in our athermal simulations. Also in figure [TBT a) 
in order to relate mobility and structure we calculate the 
SISF for different region of the material according to their 



To understand the dynamics in the soft phase and more 
generally the rheology (or mechanical response) of the 
material one would like to understand what first triggers 
the nucleation source points (local plastic events) at some 
specific locations and second how these local rearrange- 
ments interact in a cooperative manner. 

In the previous section we have analyzed the cou- 
pling between the elasticity map and the nonaffine field 
in the material and claimed that, for sufficiently large 
strain, the non-affine field (for exemple shown in fig- 
ure [T4|) is essentially irreversible (plastic) . We checked 
this assertion by comparing for each particle in the sys- 
tem the total nonaffine displacement and the purely ir- 
reversible displacement. The irreversible displacement 
field is defined as the residual displacement field result- 
ing when after each macroscopic elementary strain in- 
crement Se ^ 5.10~^ (in figure [J) the 'virtual' reverse 
shear —Se is applied on the system. For a purely re- 
versible deformation this field should cancel exactly, but 
here one observes that, for most incremental deforma- 
tion steps, a non- vanishing residual irreversible displace- 
ment field is present. With this numerical protocol we 
can extract for each particle i of the system the purely 
irreversible displacement ArJ^.^^.^, from the total non- 
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FIG. 16: On the left the self intermediate scattering func- 
tion (SISF) Fs(k,e) is plotted for k= kx — kp (squares) and 
k= ky = kp (diamonds) and for regions of different softness ci 
(thin lines, from bottom to top the rigidity ci is increased). 
Inset in the left: the relaxation strain is shown as a func- 
tion of the softness at the origin. On the right panel the 4 
point correlation function is represented for rigid (squares) 
and Lees-Edwards boundary conditions (circles). For small 
strain increments the linear behavior is highlighted. 



affine displacement ArJ^^^. Hence for each particle i the 
total nonafhne displacement over a macroscopic strain 
Ae = J2n can be decomposed into Ar^jj(Ae) = 
Arirrev('5en) + I]« ^^ev (^C") ■ The relative error 
II^Ca(Ae) - A4rev(^e)ll/II^Ca(Ae)|| averaged over 
all particle is obtained and being less than 5% confirms 
our assumption that the nonaffine displacement is dom- 
inated by an irreversible plastic contribution. Based on 
this observation we describe here the link between plas- 
ticity in the material and the local elasticity map. To 
obtain this information we study here the dynamics of 
the local rigidity ci calculated on each particle over a 
strain range of 10%. 

Figure[T7|illustrates the evolution of the coarse-grained 
modulus ci on a particle that experiences plastic activ- 
ity (i.e. that rearranges over the strain range ). W e have 
superimposed the quantity Dbf defined in [20|, which 
evaluates the degree of local deviation from affinity. This 











£ 




^sliift 


span 

1 




• \| ikl 








* lL_(u... 








III 



0.005 0.01 0.015 0.02 0.025 



e 




2,5 5 7,5 !0 12,5 15 17,5 20 22,5 25 



FIG. 17: (a): Evolution of Ci of a particle in a plastic zone. 
The dotted line represents the parameter Dbf that marks the 
local deviation from affinity. esMft and Sapan are defined for 
the point marked in red. (b): cumulative plastic activity for 
~ 0.5% and for diflferent values of tshift- 



parameter was found to very accurately distinguish be- 
tween plastic like zone {Dbf 1) and normal elastic 
zones {Dbf < 10^^). Figure [TSfa) shows the average 
behavior of the same quantity ci , averaged for all plastic 
events. The result which is typical of the dynamics of 
ci shows that before a plastic event occurs on the site 
the modulus ci decreases over a typical strain interval 
of about 0.2% to become zero or even negative at the 
plastic irreversible event, where the non-affine displace- 
ment becomes important {Dbf ^ !)• Then the local 
structure is relaxed and the local modulus gets a higher 
value (ci ~ 18 after the event). Note that this value 
ci ~ 18 is smaller than the macroscopic value for 2fj, but 
corresponds to the average value ci of the shear modulus 
at the scale W = 5 oi description. Just before a plas- 
tic event, the local shear modulus is very low ci << cT. 
The average decrease of ci before a plastic event is fit- 
ted approximately by an exponential decay to its limit 
value. We have no explanation for the moment for this 
exponential fit. It shows a characteristic strain ^ 0.2%. 

It is interesting to compare in our system this order 
parameter with other possible predictors of plasticity 
introduced in the literature. These parameters are the 
local stress [1^ [2^, HO], the local deformation strain 
[H, HI] , the local free volume [lO, HH, [111 , coordinance 
defects [s^ as well as other local criteria derived 
from macroscopic mechanics such as Tresca local yield 
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FIG. 18; Evolution of different fields around a plastic event 
and averaged over all plastic events, (a) local shear modulus, 
(b) local shear stress, (c) local density and (d) number of 
Voronoi neighbors. 



criterion or a local Mohr-Coulomb [8J, |85|- In figure 
1181 we have summarized the evolution of some of these 
fields measured locally in our model glass at the sites 
that undergo plastic rearrangement before, during and 
after the relaxation takes place. The curves are averaged 
over all plastic events over the strain studied 1000 
events). It shovifs that the density and the inverse 
compressibility are not affected by the plastic event. 
The shear stress is affected since it decreases suddenly 
after the plastic event occurred; but before the plastic 
event, the variation of shear stress is very small, and on a 
very small strain-range. Thus the conclusion appearing 
from the above analysis is that the local shear modulus 
ci is the best criterion in our Lennard- Jones glasses, to 
identify zones that are about to rearrange. 

In order to identify the effect of a very low initial local 
shear modulus ci << ci on the plastic activity over a 
larger strain span, we have calculated the shifted plastic 
activity A{i, emacro] ^shift, <^span) measured at each site i 
and every step in the macroscopic strain Cmacro and de- 
fined as the integral over a strain range espan of typically 
a few % of the local quantity Dgp starting from a shifted 
value of the macroscopic strain emacro+ishi ft from the ac- 
tual macroscopic total strain where ci is measured. Note 
that eshift and espan are adjustable parameters. 

In figure [T7] we illustrate the meaning of this two pa- 
rameters, espan is the range over which the plastic ac- 
tivity is recorded, in the limiting case of egpan — *■ oo 
one should obtain the average activity of the glass for- 
mer, eshift allows to substract a systematic bias associ- 
ated with the conditional probability to have an increased 
plastic activity for low initial values of ci . In figure [T7] 
we also plot the correlation between the average activity 
and Ci for different espan- The plastic activity is a num- 
ber that counts the number of significant plastic events. 
It is incremented by 1, as soon as Dbf > 0.2. It ap- 
pears that the probability to encounter a plastic event 
is larger for originally soft regions (ci < cT). Moreover, 
the predictive character of the structural softness of the 
material on the subsequent plastic activity of the glass 
former holds even for relatively large eshift ~ 4%. 



V. CONCLUSION AND PERSPECTIVES 

We have shown in this paper that the elastic response 
of a 2D Lennard- Jones glass at very low temperature is 
heterogeneous. We have characterized the local elastic re- 
sponse by an extensive study of the local elastic moduli 
as a function of the coarse-graining scale. We have shown 
that at a very small scale {W < 5) the elastic response 
deviates from Hookes's law, and that on an intermediate 
scale (5 < < 20) the system is not homogeneous and 
not isotropic, but becomes approximately homogeneous 
beyond this scale. However, it is difficult to identify a 
well defined characteristic scale, because all the quanti- 
ties probed in this case are power-law dependent of the 



16 



coarse-graining scale. By considering the spatial oscilla- 
tions visible in the autocorrelation function of the local 
shear modulus for W < 10, we have decided to compute 
the local moduli at the scale W = 5 where Hooke's law 
applies, and where the heterogeneities are clearly identi- 
fied. 

We have shown that these heterogeneities can be re- 
lated with the local dynamics of the particles. Regions 
with very low local shear modulus give rise to enhanced 
non-affine displacements and plastic activity. Moreover, 
the relaxation of low moduli during plastic deformation of 
the sample is faster that the relaxation of rigid zones, and 
the system seems to be decomposed into fast evolving re- 
gions with low local shear modulus (ci < cT), and rigid, 
frozen regions with high local shear moduli (ci > cT). 
The occurrence of large plastic events seems to be as- 
sociated with a large proportion of soft regions in the 
system (figurjlTj). 

We have also shown that the local shear moduli can 
be used as a predictive local criterion for plasticity. This 
is true on very short strain scales, as we have seen in 
that the lowest local shear modulus evolves in a system- 
atic way before a plastic event, with a decrease taking 
place over a typical strain of e ~ 0.2% before a plas- 
tic event. Interestingly, this is also true on larger strain 
scales, as a site with an initially low modulus tends to 
remain plastically active over long periods of strain, of 
the order of a few percent. However it is more difficult 
to identify the type of plastic events in which a site will 
be involved, which may vary considerably in size from 
isolated quadrupolar events to elementary shear bands 
spanning the whole system. The present work is a first 
step in order to predict plastic activity, and to relate it 
with the heterogeneous elastic structure in the system. 
Further studies will be needed, in order to make contact 
with existing models of plasticity and rheology of amor- 
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